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Abstract. In this paper, we present applications of the calculus developed in [9], and ob- 
tain an exact formula for the moments of random quantum channels whose input is a pure 
state thanks to gaussianization methods. Our main application is an in-depth study of the 
random matrix model introduced by Hayden and Winter and used recently by Brandao, 
Horodecki, Fukuda and King to refine the Hastings counterexample to the additivity conjec- 
ture in Quantum Information Theory. This model is exotic from the point of view of random 
matrix theory, as its eigenvalues obey to two different scalings simultaneously. We study its 
asymptotic behavior and obtain an asymptotic expansion for its von Neumann entropy. 



1. Introduction 

In the paper [9] we developed a calculus allowing to compute any moments of random 
quantum channels. It already proved useful to understand the random matrix models involved 
in the additivity violation theorems and to give improvements on lower bounds of dimensions 
needed to obtain violation of the additivity of entropy estimates (developped in [9l[T0]). In 
the present work, we study two more applications of our calculus, to new random matrix 
models introduced for quantum information theoretic purposes. 

The first application is of theoretical interest and of non asymptotic nature: we extend our 
calculus to Gaussian matrices and show that it yields explicit formulas for the moments of 
Wishart matrices and of outputs of quantum random channels. The formulae are of purely 
combinatorial nature and allow to bypass Weingarten calculus, whose asymptotic estimates 
can be involved. For this we use a 'gaussianization' method. 

The second application is a study at length of the random matrix model that was introduced 
by Hayden and Winter in [19] and used recently in [13|. IT^ Uj to refine the results of Hastings 
|18j . As a motivation, let us recall the quantum information theoretic context of this random 
matrix. A quantum channel is a linear completely positive trace preserving map ^> from 
AiniC) A4fc(C). A density matrix is a selfadjoint positive matrix of trace 1. Let = 
{x G I X]i=i ~ be the {k — l)-dimensional probability simplex. The Shannon entropy 
of X is defined to be 

k 

H{x) = - y^Xj log Xj. 

i=l 

These definitions are extended to density matrices by functional calculus: 



H{p) = — Tr p log p. 
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For a quantum channel ^ : A^n(C) — > A4k{C), its minimum output entropy is defined by 

= min Hmp)). 

peA4„{C) 

p>0,Trp=l 

The additivity conjecture for minimum output entropies was arguably one of the most im- 
portant in quantum information theory, and it can be stated as follows: 

Conjecture 1.1. For all quantum channels <I>i and ^2, one has 

(1) -f^minl^'l «) $2) = i^mm(^l) + i^min($2). 

This conjecture was disproved by Hastings in [18] as follows: 

Theorem 1.2. There exists a counterexample to the conjecture for the choice = ^>2- 

In the proof of [18j, one reason why $1 = $2 yields a counterexample is that it ensures 
that the largest eigenvalue of outputs of well-chosen inputs -Bell states- is much bigger than 
the others eigenvalues. The counterexamples to the additivity conjecture obtained so far use 
a random matrix models which we redefine in Section \6.3\ and call Zn- The main result of 
this paper is as follows: 

Theorem 1.3. The eigenvalues Ai ^ • • • ^ A„2 of Zn satisfy: 

• In probability, cnXi 1. 

• Almost surely, ^^^^ ^^^=2 ^c^n^Xi converges to a Marchenko-Pastur distribution of pa- 
rameter (? . 

• Almost surely, 

^ |21og(cn)-f + 0(1) if 0<c<l, 

as n ^ 00, where H is the von Neumann entropy. 

The interest of this result is that it yields improvements to the results of [131 [T^ [H |T8^ , as 
the only data that these papers were using was a lower bound on the largest eigenvalue of Z„, 
whereas the above theorem gives a full understanding of the eigenvalue behavior of Zn ■ 

In addition, the matrix model Zn has the novel particularity that it has two different regimes 
for its eigenvalues (one in and one in n~^). As far as we know, it is the first model in 
random matrix theory whose eigenvalues have two regimes simultaneously. 

The proof of the main theorem uses a mix of moment methods and functional calculus 
methods. It is very instructive, as the moment method is used to prove the convergence in 
distribution of the eigenvalues of smaller decay, and this goes beyond the standard intuition 
that moment methods rather give results about the larger eigenvalues. Actually, our The- 
orem 16.101 shows new kind of cancellation properties, going beyond those which are usually 
expectable with standard 'moments-cumulants' and 'connectedness' arguments. 

This paper is organized as follows. We first recall known things about Wick calculus, 
Weingarten calculus and non commutative and free probability theory. We also recall our 
graphical calculus introduced in [9] and extend it to Gaussian graphical calculus. We use it 
to obtain new non- asymptotic results for the moments of some single random channels. We 
obtain further asymptotic results in the single random channel setting, and then we come back 
to the random matrix model introduced in the bi-channel setting by Hayden and Winter, and 
compute the asymptotics of the subleading eigenvalues. 
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2. Wick calculus and Weingarten calculus 

In this section we recall known results allowing to compute expectations against Gauss- 
ian measures and Haar measures on unitary groups, as well as some standard facts in free 
probability theory. 

2.1. Wick calculus. A Gaussian space y is a real vector space of random variables with 
moments of all orders, such that each of these random variables are centered Gaussian distribu- 
tions. Such a Gaussian space comes up with a positive symmetric bilinear form (x, y) E,[xy]. 
Gaussian spaces are in one-to-one correspondence with euclidean spaces, and that isomor- 
phism of Gaussian spaces corresponds to the notion of isomorphism of euclidean spaces. In 
particular, the euclidean norm of a random variable determines it fully (via its variance) and 
if two random variables are given, their joint distribution is determined by their angle. The 
following is usually called the Wick Lemma: 

Lemma 2.1. Let V be a Gaussian space and xi, . . . ,Xk be elements inV . If k = 21 + 1 then 
E[xi ■ ■ ■ Xk] = and if k = 21 then 

I 

(2) E[xi---x,)]= j2 n^^[^^™^^-j 

P={{h,ji},-,{ii,ji}} "1=1 
pairing of {l,...,k} 

In particular it follows that if xi, . . . ,xP are independent standard Gaussian random variables, 
then 

p 

E[x^...xJ''] =n(2fc.)!! • 

1=1 

For a proof, see for instance [27]. It is possible to extend the notion of Gaussian space to a 
complex Gaussian space. A complex valued vector space V is called a Gaussian space if and 
only if, for any real structure on V, the pair (Re{V),lm{V)) is a real valued Gaussian space. 
One checks readily that in the case of a complex Gaussian space, the Wick lemma 12.11 holds 
with the exact same statement. 

We will usually denote by Gn,m (or G when there is no ambiguity) the standard complex 
Gaussian random matrix n x m. It has the distribution exp(—NTr{GG*))dG where dG is 

the Lebesgue measure on the n x m complex matrices properly rescaled, and G* = G is the 
standard operator algebraic adjoint. 

Since in this paper we shall mostly be concerned with traces of products of random matrices, 
we need to introduce one last notation for generalized traces which we borrow from [6]. 
For some matrices Ai,A2,...,As £ A4„(C), some permutation a £ Sp and some function 
t : {1, . . . , p} —f {1, . . . , s} we define 

Tr.,(Ai,...,^)= n Tr(n 

When s = p, we use the simplified notation Tro-,t(Ai, . . . , Ap) = Tro-^id(Ai, . . . , Ap). We also 
putTr^iA) = Ti^{A,A,...,A). 

2.2. Weingarten calculus. In this section, we recall a few facts about Weingarten calculus. 

Definition 2.2. The unitary Weingarten function Wg(n, o") is a function of a dimension 
parameter n and of a permutation a in the symmetric group Sp. It is the inverse of the 
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function a ^ rt^" under the convolution for the symmetric group (^a denotes the number of 
cycles of the permutation a). 

Notice that the function a i— >■ n"^"^ is invertible as n is large, as it behaves like n^6e as 
n ^ oo. We refer to [12] for historical references and further details. We shall use the 
shorthand notation Wg(c7) = Wg(n, a) when the dimension parameter n is obvious. 

The function Wg is used to compute integrals with respect to the Haar measure on the 
unitary group. 

Theorem 2.3. Let n be a positive integer and {ii, . . . ,ip), {i[, . . . , i'p), (ji, . . . , jp), . . . ,jp) 
be p-tuples of positive integers from {1,2,... , n}. Then 



(3) / U,,,,---U,^,^U,^^>---U,^,^dU 

JU{n) 



E -^nC,,) • • • ■ ■ ■ Kk,,, Wg(n, ra-i). 



If p ^ p' then 



(4) / Ui,,, 

Juin) 



_ ■■■U^,3^^,,[■■■U,,,>, dU = 0. 
U{n) ' ' f f 

We are interested in the values of the Weingarten function in the limit n — > oo. The 
following result encloses all the data we need for our computations about the asymptotics of 
the Wg function; see [8] for a proof. 

Theorem 2.4. For a permutation a G Sp, let Cycles(cj) denote the set of cycles of a. Then 

(5) Wg(n,a) = (-ir-#'^ W Wg(n,c)(l + 0(n-2)) 

cSCycIes(o-) 

and 

(6) Wg(n,(l,...,d)) = (-l)'^-ic,_i W (n-j)-' 

where Ci = ^^^^|| ^| is the i-th Catalan number. 

A shorthand for this theorem is the introduction of a function Mob on the symmetric group, 
invariant under conjugation and multiplicative over the cycles, satisfying for any permutation 
a G Sp\ 

(7) Wg(n, a) = n-(P+l'^l)(Mob((j) + 0(n-2)). 

where \a\ = p — i^cr is the length of a, i.e. the minimal number of transpositions that multiply 
to a. We refer to 22] for details about the function Mob. 



2.3. Elementary reminder of non-commutative and free probabiUty theory. A non- 
commutative probability space is an algebra A with unit endowed with a tracial state ip. An 
element of A is called a (non-commutative) random variable. In this paper we shall be mostly 
concerned with the non-commutative probability space of random matrices (A1„(L°^"(ri,P)),E[n-i Tr(-)]) 
(we use the standard notation L°°-{n,F) = np^iLP{n,F)). 
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Let (ai, . . . , afc) be a A; -tuple of selfadjoint random variables and let C{Xi,...,Xk) be 
the free *-algebra of non commutative polynomials on C generated by the k indeterminates 
Xi, . . . , Xk- The joint distribution of the family {ai}^^i is the linear form 

/^{ai,...,afc) : . . . ,Xk) C 

P ^ v?(P(ai, . . . ,afc)). 

Given a /c-tuple (ai, . . . , a^) of free random variables such that the distribution of a, is /Xq. , 
the joint distribution /^(ai,...,afe) is uniquely determined by the //a/s. A family (a", . . . , a'^)n of 
/c-tuples of random variables is said to converge in distribution towards (oi, . . . , a^) iff for all 
P G C{Xi, . . .,Xk), lJ'{a^,...,ai){P) converges towards ^(ai,...,afc)(-P) as n ^ oo. 

The following result is contained in [23j and is of crucial use for us. 

Lemma 2.5. The function d{a,T) = \a~^T\ is an integer valued distance on Sp. Besides, it 
has the following properties: 

• the diameter of Sp is p — 1 ; 

• d{-,-) is left and right translation invariant; 

• for three permutations (Ti,(T2,t G Sp, the quantity d{T,ai) + ci(r, (T2) has the same 
parity as d{ai,a2); 

• the set of geodesic points between the identity permutation id and some permutation 
(T ^ Sp is in bijection with the set of non-crossing partitions smaller than tt, where 
the partition tt encodes the cycle structure of a . Moreover, the preceding bisection 
preserves the lattice structure. 

We finish by gathering the bare minimum of free probability theory needed towards the 
main results of this paper. We skip the definition of freeness, as we won't need it. Free 
cumulants are multilinear maps indexed by non-crossing partitions a G NC (p) on p elements 

: ^4 X • • • X ^4 ^ C 

p times 

such that 

(8) ^ K^{xi, . . . ,Xp) =E^[xi, . . . ,Xp] 

ni^a£NC{p) 

for all non-crossing partitions a G NC (p) , where Eq- [ the product over the blocks 

{xij^, . . . ,Xi.} of fj, oi E,[xi-^ . . . Xi-]. Cumulants are known to be multiplicative over blocks, 
therefore a special role is played by the cumulant corresponding to the maximal partition Ip, 
which we denote by K{ai, . . . , ap) := Kip(ai, . . . , ap). 

We will need free cumulants for computational purposes, in order to identify free Poisson 
distributions. Let us mention for the interested reader that the main property of the free 
cumulants is that mixed cumulants of free variables vanish. 

We recall that the free Poisson distribution of parameter c is given by 



TTc = max(l - c,0)So + — Mi+c-2V^,i+c+2^{x) dx. 

It is characterized by the fact that all its free cumulants are equal to c. Although we will 
not need this fact, it is worth to mention that it has a semigroup structure with respect 
to the additive free convolution of Voiculescu (see, e.g. 123]). It is also sometimes called 
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Marchenko-Pastur distribution. One can compute (minus) the entropy of this probability 
distribution 



(9) Kc = J xlog X diTcix) 



^ + clog c if c ^ 1, 
# if < c < 1. 



3. Unitary and Gaussian graphical calculi 

In this section we recall briefly the results of [9| for the convenience of the reader and in 
order to make the paper self contained. Then we introduce the Gaussian graphical calculus 
and we present a first application of it to Wishart matrices. 

3.1. Axioms of unitary graphical calculus. The purpose of the graphical calculus intro- 
duced in [9] is to yield an effective method to evaluate the expectation of random tensors 
with respect to the Haar measure on a unitary group. The tensors under consideration can 
be constructed from a few elementary tensors such as the Bell state, fixed kets and bras, and 
random unitary matrices. In graphical language, a tensor corresponds to a box, and an ap- 
propriate Hilbertian structure yields a correspondence between boxes and tensors. However, 
the calculus yielding expectations only relies on diagrammatic operations. 

Each box B is represented as a rectangle with decorations on its boundary. The decorations 



are either white or black, and belong to S{B) U S* (B). Figure 1(a) depicts an example of box. 



M * 



(b) 



M p-^ y 

(c) 



c A 


• 


• 






c B 


■ □ 


I 



(d) 



Figure 1. Basic diagrams and axioms 



It is possible to construct new boxes out of old ones by formal algebraic operations such as 
sums or products. We call diagram a picture consisting in boxes and wires according to the 
following rule: a wire may link a white decoration in S{B) to its black counterpart in S*{B). 
A diagram can be turned into a box by choosing an orientation and a starting point. 

Regarding the Hilbertian structure, wires correspond to tensor contractions. There exists 
an involution for boxes and diagrams. It is antilinear and it turns a decoration in S{B) into its 
counterpart in S*{B). Our conventions are close to those of [3 [20], and are hoped to be familar 
to the reader acquainted with existing graphical calculus of various types (planar algebra 
theory, Feynman diagrams theory, traced category theory). Our notations are designed to fit 



well to the problem of computing expectations, as shown in the next section. In Figure 1(b) 



1(c) and 1(d) we depict the trace of a matrix, multiplication of tensors and the tensor product 



operation. For details, we refer to [9]. 
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3.2. Planar expansion. In this subsection we describe the main appUcation of our calculus. 
For this, we need a concept of removal of boxes U and U. A removal r is a way to pair 
decorations of the U and U boxes appearing in a diagram. It therefore consists in a pairing 
a of the white decorations of U boxes with the white decorations of U boxes, together with 
a pairing /3 between the black decorations of U boxes and the black decorations of U boxes. 
Assuming that D contains p boxes of type U and that the boxes U (resp. U) are labeled from 
1 to p, then r = (a, /?) where a, (3 are permutations of Sp. 

Given a removal r G Rem(I?), we construct a new diagram associated to r, which has the 
important property that it no longer contains boxes of type U or U . One starts by erasing the 
boxes U and U but keeps the decorations attached to them. Assuming that one has labeled 
the erased boxes U and U with integers from {1, . . . one connects all the (inner parts of 
the) white decorations of the z-th erased U box with the corresponding (inner parts of the) 
white decorations of the a(i)-th erased U box. In a similar manner, one uses the permutation 
/3 to connect black decorations. 

In [9], we proved the following result: 

Theorem 3.1. The following holds true: 

3.3. Gaussian planar expansion. Now, we consider the case where in our diagrams we 
allow a new special box G corresponding to a Gaussian random matrix. We shall address the 
same issue as in the unitary case: computing the expected value of a random diagram with 
respect to the Gaussian probability measure. 

To start, consider D a diagram which contains, amongst other constant tensors, blocks 
corresponding to independent Gaussian random matrices of covariance one (identity). One 
can deal with more general Gaussian matrices by multiplying the standard ones with constant 
matrices. Note that a block can appear several times, adjoints of blocks are allowed and the 
diagram may be disconnected. Also, Gaussian matrices need not to be square. 

The expectation value of such a random diagram D can be computed by a removal pro- 
cedure as in the unitary case. Without loss of generality, we assume that we do not have in 
our diagram adjoints of Gaussian matrices, but instead their complex conjugate blocks. This 
assumption allows for a more straightforward use of the Wick Lemma 12.11 As in the unitary 
case, we can assume that V contains only one type of random Gaussian blocks G; the other 
independent random Gaussian matrices are assumed constant at this stage as they shall be 
removed in the same manner afterwards. 

A removal of the diagram P is a pairing between Gaussian blocks G and their conjugates 
G. The set of removals is denoted by RemG'(D) and it may be empty: if the number of 
G blocks is different from the number of G blocks, then RemG(X>) = (this is consistent 
with the first case of the Wick formula ([2|)). Otherwise, a removal r can identified with a 
permutation a G Sp, where p is the number of G and G blocks. Let us stress here the main 
difference between the notion of a removal in the Gaussian and the Haar unitary case. In the 
Haar unitary (or the Weingarten) case, a removal was associated with a pair of permutations: 
one had to pair white decorations of U and U boxes and, independently, black decorations of 
conjugate blocks. On the other hand, in the Gaussian/Wick case, one pairs conjugate blocks: 
white and black decorations are paired in an identical manner, hence only one permutation 
is needed to encode the removal. 
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paired boxes 



Figure 2. Pairing of boxes in the Gaussian case 
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Figure 3. Diagram of a Wishart matrix 



To each removal r associated to a permutation a G iSp corresponds a removed diagram V^. 
constructed as follows. One starts by erasing the boxes G and G, but keeps the decorations 
attached to these boxes. Then, the decorations (white and black) of the i-th G block are 
paired with the decorations of the a(i)-th G block in a coherent manner, see Figure [2j 

The graphical reformulation of the Wick Lemma [2. II becomes the following theorem, which 
we state without proof. 

Theorem 3.2. The following holds true: 

3.4. Moments of Wishart matrices. As a first application of our Gaussian graphical cal- 
culus, we compute the moments of traces of products of Wishart matrices. By definition, a 
Wishart matrix of parameters (n, h) is a positive random matrix W E A^„(C) such that 

W = G-G\ 

where G € A^nxfc(C) is a standard Gaussian random matrix. In our graphical formalism, 
since we only consider Gaussian random matrices, the previous equation corresponds to the 
graphical substitution in Figure [3l round decorations correspond to n-dimensional complex 
Hilbert spaces C" and square-shaped labels correspond to C*^. 

The same problem of computing expected values of traces of Wishart matrices was consid- 
ered in [6l [T5l [T71 [2T] and we shall re-derive Corollary 3 of Theorem 2 from [6] . The general 
covariance case (Theorem 2 in [6]) can be easily derived from the result below. 

Proposition 3.3. Let Wi,W2, ■ ■ ■ ,Ws he independent Wishart matrices of unit covariance 
and parameters {n,ki), {n,k2), {n,ks) respectively. For a permutation a & Sp and a 
function t : {1, . . . , p} ^ {1, . . . , s} , the following holds true: 

s 

(10) E[TT,,t{wi,...,ws)]= n^r^^*^''''"^' 

aeSp{t)j=l 

where Sp{t) = {a Sn\t = t o a}. Every permutation a G <Sp{t) leaves the level sets of t 
invariant, and it induces on each set t~^(j) a permutation aj (j = 1, . . . ,s). 
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Figure 4. Monomials of traces of Wishart matrices 

Proof. We consider the diagram T) corresponding to the left hand side of equation (jlOh . It 
contains n Wishart boxes from the set {VFi, . . . , Wg} which are wired according to the per- 
mutation a (see Figure [4]). Computing the expectation of the diagram T) is rather straightfor- 
ward using our graphical calculus. Since we are dealing with s independent Gaussian matrices 
Gi, . . . ,Gs (recall that Wj = GjG*) one needs to apply Theorem 13.21 s times, once for each 
Gaussian matrix Gj. Each box Gj appears |t~^(j)| number of times and, using Theorem 13.21 
we get 

where each permutation aj G encodes the removal procedure for the Gj boxes. 

Diagrams obtained after the successive removal procedures I'oi,...,^^ are made of loops of 
two types: loops associated to the n-dimensional space C" and loops associated to "internal 
spaces" C^^ . In order to count the number of loops of each dimensionality, let us first notice 
that the set of s-tuples of permutations (ai, . . . , a^) is in bijection with the set of permutations 
a G Sp{t) defined in the statement of the theorem. 

For such a permutation a G Sp(t), let us count the number of loops corresponding to traces 
over C'^^ . Initially, the pj decorations of the Gj boxes are connected in the simplest manner: 
the kj decoration of the i-th Gj box is connected to the corresponding decoration of the Gj 
box with the same index i. The j-th removal procedure, encoded by the permutation aj 
produces then a number of #(id~^ aj) = ^Uj loops. Hence, the contribution of the C'^J-type 
loops is k- . 

The computation of the loops associated with C" is more involved, since the decorations 
are already non-trivially linked by the permutation a. Since a may not respect the level sets 
of the function t, one needs to consider the global action of a, the restrictions Uj not being 
sufficient in this case. Since the boxes are initially connected by a and the removal procedures 
add wires according to the permutation a, the total number of loops is ^{a~^a). Adding all 
loop contributions, one obtains the announced formula (jlOp . □ 



Remark 3.4. One can consider in the graphical model more general covariances and obtain 
Theorem 2 of ^ in its full generality. All there is to be done is to add constant tensors 
associated with covariance matrices in our diagrams. After the successive removal procedures, 
one is left with loops and traces of monomials in these constant matrices. Since our purpose 
in this section was to illustrate the Gaussian graphical calculus, we leave the details of this 
more technical generalization to the interested reader. 

4. Application of Gaussianization: pure states through random quantum 

channels 

4.1. Random single channel model. In this section we present an important application 
of the Gaussian diagrammatic calculus: we compute eigenvalue statistics for the action of a 
random quantum channel on a pure quantum state. By definition, a quantum channel $ : 
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U 



Y 



cX 



U* 



Figure 5. Diagram for a quantum channel 



<57 



Figure 6. An equivalent diagram for quantum channels with rank one X and Y 



M.n{C) 7V4„(C) is a trace preserving, completely positive map. According to Stinepsring 
theorem, such a linear application can be written as 

= = Trfc [U{X (g) Y)U*] , 

where C/ is a unitary matrix in U{nk) and y is a /c-dimensional rank-one projector. A 
diagrammatic picture of the above formula is presented in Figure [5l The set of quantum 
channels can be endowed with a natural probability measure by fixing the projection Y and 
picking U uniformly, with respect to the Haar measure on the unitary group U{nk). This 
is the model of randomness we refer to when we speak or random quantum channels and it 
has received a lot of attention from the quantum information community [9l |T9] . From the 
definition of one can see that the Weingarten calculus developed in ^ may be applied 
to this situation, since random unitary matrices are a key-element in the problem. However, 
when random quantum channels are presented with rank-one inputs (or pure states) , we show 
that one can use the simpler Gaussian calculus. Using this approach, we shall recover some 
exact formulas for the moments of the output from [9j, as well as some asymptotic results 
from [22j . 

We are interested in the output random matrix 



(11) 



where X is a rank-one projector. The main result, obtained in [22], is as follows: 



Proposition 4.1. Let W 



G-G* & ^AniC) be a Wishart matrix of parameters {n,k). Then 
Z = $(X) = W/Tt{W). 



Observe that this result does not depend on the choice of X, Y due to the invariance of the 
Haar measure. 

The main point is that one can show (see [22]) that the eigenvalues of Z, that is the 
normalized eigenvalues of W, are independent on the trace of W. This implies that we we 
can apply the results on Wishart matrices developed in Section [3.41 to this particular case. 



4.2. Exact moments. We provide in this section exact formulas for the moments E[Tr(ZP)] 
of the output of a random quantum channel. Other formulas for the same quantities (as well 
as some recursion relations) have been obtained in [22^ [25l I26j . 
Using the Gaussianization trick, we have that 

E[Tv{WP)] 



E[TliZP)] 



E[Tv{W)P] ' 
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where is a Wishart matrix of parameters (n,k). One uses Proposition 13.31 to compute 
E[Tr(M^P)] and E[Tr(VF)P]: 



E[Tr{WP)] = k*^n'' 



E[Tr{W)P] = ^ 

aeSp 

where "y = {p p — ^ • • • 2 1) € 5p is the full cycle. In the second formula above, one recognizes 
the generating polynomial of the number of cycles of permutation of p objects evaluated at 
nk. This is known to be equal to nk{nk + l){nk + 2) ■ ■ ■ (nk+p — l), and one gets the following 
theorem. 



Theorem 4.2. 



'p-i 



(12) E[TY(ZP)] = Ylink + ^*^'^ 



j=0 / a&Sp 



This is exactly as formula (10) from [9J, which was obtained via Weingarten formula. The 
approach followed here is more straightforward and does not use unitary integration. It is 
based on the purely combinatorial Wick formula and the Gaussianization trick. 



4.3. Asymptotics. We now look at the probability distribution of the output random matrix 
Z when one (or both) of the parameters n and k grow to infinity. The asymptotic behavior 
of random matrices has been one of the main objects of study in random matrix theory; 
for instance, it is in this large dimension regime that the freeness phenomenon appears. 
In the particular case of random quantum channels under study here, this question has an 
interesting physical motivation: large dimensional Hilbert spaces model physical systems with 
a large number of degrees of freedom. This point of view has been discussed in the quantum 
information theory literature (see O [26l [22l [2] ) . Although some of what follows has already 
been treated in [22], the approach of this paper has the merit of being self-contained and 
illustrates perfectly the power and range of the Gaussian graphical calculus. 

We split the results according to three possible asymptotic regimes, depending on which 
of the parameters n and/or k is large. Of special interest is the third regime, when both 
parameters grow to infinity, but at a constant positive ratio c > 0. 

Theorem 4.3. Let Z = <I>^'-^(A) the output of a random quantum channel ^, where X and 
Y are rank-one projectors. 

(I) In the regime n fixed, k oo, the limiting spectral distribution of Z is almost surely 

(II) In the regime k fixed, n — > oo, Z tends almost surely to a variable that has eigenvalues 

1/k with multiplicity k and with multiplicity n — k. 
(Ill) In the regime n,k ^ oo, /c/n — > c > 0, cnZ converges almost surely to a free Poisson 
distribution of parameter c. 

Proof. In the first regime, 

E[Tr(ZP)] -(nA;)-P ^ A:#°n#(^"'"). 
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Permutations a which give non vanishing contributions are those such that #a = p, hence 
a = id. At the end, we obtain 

hm E[Tt{ZP)] = n^-P, 

hence the hmiting spectral distribution of Z is 5i/n- 

In order to prove the almost sure convergence, we show that the empirical measures 



1 " 



i=l 



converge almost surely to the limit (5x/„ (which is equivalent to the fact that almost surely, ev- 
ery eigenvalue of Z converges to 1/n - recall that n is fixed). As usual, almost sure convergence 
of moments suffices and we set our goal to prove that for all p 



a.s. lim Tr{ZP) = n^'P. 

fc— »oo 

A standard application of Chebyshev's inequality and Borel-Cantelli's lemma shows that it is 
enough to verify that for all integers p the series of variances is summable: 



oo 

[(Tr(ZP) -ETr(ZP))^ 



k=l 



< OO. 



Let us compute separately E[Tr(ZP)^] and E[Tr(ZP)]^ using formula ()12p . For the first expec- 
tation, one needs to introduce the permutation 



(13) 72 = {p{p- 1) 

One then has 



2 1)(2p(2p-l) ••• (p + 2) {p+l))eS2p. 



E[Tr(Z 




aeS2p 

-1 



oieS2p 



The first contribution (of order k^) in the last sum is given by a = id and it is equal to ■n?~'^P 
(recall that 72 has 2 cycles). The second order in k is given by transpositions a = (ij). In this 
case, |72'^a| = 2p — 3 if z and j belong to the same cycle of 72 and |72^^a[ = 2p — 1 otherwise. 
Hence, we obtain 



E[Tr(Zf)2 



2p(2p-l) /I 
2nk ^ U2 



n 



2-2p 



+ Z {P^n'-^P + p{p - l)n3"2p) + O ( ^ 



n2-2p + ip(p _ l)ni-2p(n2 - 1) + O (^^^ 
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Using the same ideas, E[Tr(ZP)]^ is easily computed: 





/2p-l 


= 


n 




V i=o 











2nk 




aeSp 



+ - l)n'-*(n^ _ 1) + O ("^^ 



and one concludes that K[Tic{Z'P)'^ — E[Tr(ZP)]^ = 0{k^'^) and thus the covariance series 
converges, finishing the proof. 
In the second regime. 

The non- vanishing contribution is given by a = 7, and thus 

lim E[Tr(ZP)] = k^'^ . 

n— >oo 

In other words, for large n, Z has the following eigenvalues: 

• 1/k with multiplicity k\ 

• with multiplicity n — k. 

The proof of the almost sure convergence follows the same lines as in the previous case and 
it is left to the reader. 

In the third regime, after making the substitution k = cn, the asymptotics are 

(14) E[Tr(ZP)] ~ n-^Pc-P ^ c#"n#"+#(^"'"). 

a£Sp 

Since 

(15) #a + #(7~'a) = 2p - (|a| + {j-^aD ^p + l, 

one should rescale the matrix Z by a factor of n. In fact, in order to avoid some unnecessary 
complications, we shall rescale Z by cn. We get: 

Contributing permutations are those for which we have equality in equation (|15p . that is 
\'^\ + |7~^ck| = ItI = P ~ 1- These are permutations on the geodesic id — 7 and are knonw to 
be in bijection with non-crossing partitions a E NC{p). Thus 

E[tr„((cnZ)P)] ~ ^ c#-. 

One recognizes the moment-cumulant formula from free probability theory. Hence, the limit- 
ing distribution of cnZ has cumulants of all order equal to c and one identifies the free Poisson 
distribution of parameter c. Let us now show that almost sure convergence holds: 

lim E[tr„((cnZ)^)] = d^'^ almost surely. 

aeNC{p) 
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Using the same classical technique as in the first regime, we show that the series 

Y {n{tTn{{cnzrf)]-ntTn{{cnzy)?) 

n 

converges. We start by evaluating E[(tr„((cnZ)^)'^)] up to the second order in n. Using the 
permutation 72 defined in (fT3]) and the Gaussian graphical calculus, we have 

E[(tr„((cnZ)f)2)] ~ ^ ^*a^2p-2-{\a\ + \^-^a\) ^ 

Using similar ideas as before, |q| + |72^"'^a| ^ I72I = 2p — 2, with equality iff a is on the 
geodesic between id and 72 . Given the 2-cycle structure of 72 , geodesic permutations a admit 
a decomposition a = a' + a" , where a' G 2, . . . ,p} = Sp and a" G S{p+l,p+2, . . . , 2p} ~ 
Sp are themselves geodesic permutations id ^ q' ^ 7 and id a" — > 7. Of course, in this 
case, #a = #a' + #a" and thus 

EiitT^iicnZrf)] ^ Y = ( E 

id— >o"— >7 

By a standard parity argument, the function S2p B a {\a\ + [72^ ""^al) mod 2 is constant and 
thus there is no term in the asymptotic development of E[(tr„((cnZ)^)^)]: 

E[itVn{{cnZ)Pf)]=l Y c*") +0(n-2). 

\ id^ci^7 / 

Similar ideas applied to the formula (fHI) yield the same conclusion: 

E[tini{cnZ)P)] = Y c#" + 0(n^2)_ 

Taking the square of this last equation and comparing with the previous one, we conclude 
that the general term of the covariance series behaves asymptotically as 0{n~'^). This implies 
that the series is convergent and we conclude that the almost sure convergence holds. □ 

Even though Gaussianization results are exact and do not require the detour through 
Weingarten calculus, it is not clear how to apply them when the input is not one dimensional. 
However it is natural to wonder about the asymptotics in this case as well. The calculus that 
we introduced in [9] is crucial for that and this is the object of the Section [5j 

4.4. Almost sure convergence for entropies. In this section, we improve the almost 
sure convergence of moments into the almost sure convergence of any continuous function 
with polynomial growth. Since the set of functions that it applies to is larger, this type 
of convergence is stronger than the weak convergence. We deduce corollaries for quantum 
information theory, and the techniques developed in this section will be useful towards the 
end of the paper. The technique of proof of this result is inspired from p2]. 

Theorem 4.4. Let f be a continuous function on M with polynomial growth and let be 
a sequence of probability measures which converges in moments to a compactly supported 
measure v. Then f f dv^ J f dv. 
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Proof. Let K > 1 be a constant such that the interval [—(K—1),K—1] contains the (compact) 
support of the Hmit measure i'. It follows that, for all integer power s ^ 0, 



(16) 



lim K-''' I x^''+'' dv{x) = 0. 



Moreover, since the measures Vn converge in moments to for all e > 0, there exists an r 
large enough such that for all n large enough. 



(17) 



K 



-2r 



For some fixed 5 > 0, the Weierstrass theorem produces a polynomial P such that \f{x) 
P{x)\ < 5 for all x G [—K,K]. Then, we have 



j fdVn- j fdu ^ J\f-P\dUn + J\f-P\dl^+ J P 



d{Vn 



Since the polynomial approximation holds on the support of the second term above is less 
than 5. Using the convergence in moments of the probability measures m„, the last term 
can be seen to be less than 5 for n large enough. We focus now on the first term above, 
j \ f — P\ dun. By the polynomial approximation, J \f — P\ dvn ^ S + f\x\>K \f ~ ^\ 
Since / has polynomial growth, one can find a constant g > such that \f{x) — P{x)\ ^ x^'^ 
for all |x| ^ K. Using the Chebyshev inequality on the last integral, we have for all r ^ 1 

J\x\^K Jr J 

The convergence in moments together with equations (I16p and (1170 imply that, for r and n 
large enough, the above expression can be made arbitrarily small, which allows to conclude. 

□ 

Corollary 4.5. Almost surely, in the limit n oo, the von Neumann entropy of the matrix 
Z satisfies 

logn-i + o(l) if c^l. 



H{Z) 



log(cn) - f + o(l) if < c < 1. 



Proof. Let us assume that c ^ 1, the other case being similar. We use Theorem 14.41 for the 
function x i— > x log x which is continuous and of polynomial growth on the domain M+ , and 
for the empirical spectral measures of the matrices cnZ. It follows that, almost surely when 



n 



oo, 



1 " f 

— > cnAj log(cnAi) = / tlogtdfiJt) + o(l), 



where Ai ^ ■ ■ ■ ^ A„ are the eigenvalues of Z. 

Simplifying this expression and using the value of the right-hand side integral from equation 
dll), we conclude: 



i/(Z) = -^AilogAi=log 



i=l 



™-^ + o(l). 



□ 
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A formula of Page [24] states that the mean entropy of a random density matrix Z^""-^^ € 
A^n(C) obtained by tracing out a A;-dimensional environment is given by (here, n ^ k are 
fixed) : 

nk ^ ^ 

E/7(z(«'^))= y 

^ ^ ^ j 2k 

j=k+i 

One could obtain a weaker version of Corollary 14.51 from Page's formula by letting n tend to 
infinity and using the dominated convergence theorem. 

5. ASYMPTOTICS OF A SINGLE RANDOM QUANTUM CHANNEL FOR GENERAL STATES 

5.1. The model. We are interested in single random quantum channels, and study the as- 
ymptotic behavior of the output of such channels for more general input states than rank one 
projectors. The Gaussian planar expansion can not be used in the more general cases, so we 
need Weingarten planar expansion. One may consider the following general model: 

(18) Tr^(X) ~ {n')*''u*^M^), 

where s, u G M are fixed parameters and a; is a random variable in some non-commutative 
probability space with trace ip. In this section, we will deal only two special cases of interest 
of the above formula. The first one is motivated by quantum information theory: X is a rank 
r projector. This choice corresponds to s = 0, « = r and x = r~^. The second special case we 
consider will seem natural to the reader with a free probabilistic background: X converges in 
moments to a non-commutative random variable x. To get this particular case from formula 
(lisp one has to put s = u = 1 (this amounts to taking a normalized trace in the left hand 
side). Note however such an input matrix is not normalized, and one has to take into account 
the trace one restriction for quantum state. 

Let us recall here the formula for the moments of the output Z = ^{X) of a random 
quantum channel (see [9]): 

(19) E[Tr(ZP)] = J2 ^*°n*^^~'"^Tr^(^)Wg(a;3-i), 
where 7 is the full cycle permutation ^ = {p p — 1 •■■ 2 1) £ Sp. 

5.2. Rank r projectors. Plugging, for all /3 G Sp, Ti[^(X) = r'^l^r^P = j--!/^! in the previous 
equation, we obtain 

(20) E[Tt{ZP)]= A;#"n#(^^'"V-I^IWg(a/3-^). 

We study, as usual, the three asymptotic regimes n fixed, k ^ 00, k fixed, n — > 00 and 
n,k ^ 00, k/n ^ c. 

Proposition 5.1. Depending on the asymptotical regime, the almost sure behavior of Z is 
given by: 

(I) When n is fixed and k ^ 00, the output density matrix Z converges almost surely to 
the maximally mixed state 

p* — 

n 

(II) When k is fixed and n — > cxo, the output density matrix Z, restricted to its support of 
dimension rk converges to l/{rk)lrk; 
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(III) Finally, in the third regime k/n c, the empirical spectral distribution of the matrix 
rkZ converges to a free Poisson distribution of parameter rc. 

Proof. Using the Weingarten asymptotic 'Wg{a/3~^) ~ {nk)~P~^°'^ ^1, the exponent of k in 
equation (j20p is given by ^a—p — |a/3~^|. This reaches its maximum of zero when a = /3 = id. 
Hence, to the first order in k, we have 

E[Tr{ZP)]=n^-P + o{l), 

and the conclusion follows. 

The second regime is very similar, and one gets at the end (this time up to the first order 
in n) 

E[Tr(ZP)] = {rk)^~P + o{l). 

As for the third regime, making the substitution k = cn, we obtain the following asymptotic 
relation: 

E[Tr(ZP)] - r-l^lc-(l"l+l"^"'l)n-(l"l+l^"'"l+2|"/3-^l)Mob(a/3-i). 

The exponent of the large parameter n in the last formula is minimized when id — > a = /? — > 7 
is a geodesic in Sp. Hence, 

E[TV(ZP)] ~ n^-P (rc)-l^lMob(a/3"i). 

id— >«— >7 

Thus, the normalized trace of the p-th power of the matrix rkZ converges to 

id^a^7 (7eNC(p) 

and one recognizes easily the moment-cumulant formula for the Marchenko-Pastur distribution 
of parameter rc (see the reminders of Section 12. 3p . 

The above results have been proved to hold for the convergence in moments. Borel-Cantelli 
techniques (see [9] for a sample) can be easily used to show that the stronger almost sure 
convergence holds in all three cases. □ 

5.3. Normalized macroscopic inputs. We now consider matrices X which have a macro- 
scopic scaling Ti{XP) ~ n ■ ip^x^), where x is some non-commutative random variable. One 
has of course to normalize such input matrices, and we shall consider 

With this normalization, the moments of the output matrix Z = ^(X) are given by 

E[Tr($(X)P)] 



E[Tr(ZP)] = E[Tr(^>(X)P)] = E 



(Tr X)P 



(TrX)P 



As in the previous section, we consider different asymptotic regimes for the integer param- 
eters n and k. However, it turns out that the regime k fixed, n ^ 00 is more involved, and 
its understanding requires some more advanced free probabilistic tools. To an integer k and 
a probability measure /i, we associate the measure defined by 

m = (1 - 1) '^0 + 

Proposition 5.2. The almost sure behavior of the output matrix Z = ^{X) is given by: 
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(I) When n is fixed and k ^ oo, Z converges almost surely to the maximally mixed state 

P* — In • 

n 

(II) When k is fixed and n oo, the empirical spectral distribution of fiknZ converges to 
the probability measure v = [/X(fc)]^'^^, where ffl denotes the free additive convolution 
operation, p. is the probability distribution of x with respect to ip: ip{xP) = J t'^ dp{t) 
and Jl is the mean of p, p = fix). 
(Ill) When n, fc — > oo and k/n — > c, the empirical spectral distribution of the matrix nZ 
converges to the Dirac mass 6i. 

Proof. We start with the simplest asymptotic regime, n fixed and k ^ oo. Plugging in the 
scaling for Tr^(X) in formula p^ . we get 

E[Tr(ZP)] ~n~P^(x)-P ^ fc#"n#(^"'")n#^(/p^(x)(nA;)-P-l"^"'l Mob(a/3-^). 

In order to find the leading term in the preceding sum, one has to minimize the exponent of 
k, \a\ + |a/5~^|. This expression attains its minimum at a = /3 = id. At the end, one finds 
E[Tr(Z^)] ~ n^~'P and concludes that the output matrix Z converges to the maximally mixed 
state p* = I„ /n. 

Let us look now at the second regime, k fixed and n ^ oo. The asymptotic moments of Z 
are given by 

E[Tr(ZP)] ^p{x)-P A;-(l"l+l"^"'l)n-(l^l+l"'^''l+l'^''"IV/3(a;)Mob(a/3-^). 

The dominating terms in the preceding sum are given by permutations such that + |a/?~^ | + 
|7~^a| is minimal. Permutations {a, 13) which saturate the triangle inequality + |a/3~^| + 
|7~^a| ^ I7I = p—1 are elements of the geodesic id — > /3 — > a — > 7 and can be put in bijection 
with non-crossing partitions cr ^ r G NC{p) using Lemma 12.51 We obtain 

-EmUpknZy)] V k'^*^-*''ip^(x)Moh(a,T). 

Using the fact that k~'^'^ipa-{x) = tpa{p{k)) ^iid applying the moment-cumulant formula (|23j. 
pp. 175), we get 

^E[TT{{pknZr)] ^ k^*^ E MPik))Moh{a,T)= ^ k^*- ^^p^^,^), 

reNCip) aeNC{p) TeNC{p) 

where k denotes the free cumulant. We conclude that the random matrix pknZ converges in 
distribution to a probability measure v which has free cumulants K,p{v) = k'^Kp{p(^j^-^), and the 
conclusion follows. 

We turn now to the third regime, where both n and k grow to infinity at a constant ratio 
c > 0. After making the substitution k = cn, one obtains the following equivalent: 

E[Tr(ZP)] ~ (^(x)-P Yl ^~^'"'+'^"'"'+'^'+''"^"''^c-(l"l+l"^"'IV/3(a;)Mob(a/3-^). 

The expression to minimize in this case is \a\ + |7~^a| + + 2\af3~^\. By the triangle 
inequality, (cf Lemma [2. 5p . the sum of the first two terms is at least [7! = p — 1 and the other 
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terms are positive; hence, the (negative) exponent of n is at least p — 1, and the bound is 
reached for a = /3 = id. To the first order in n, the asymptotic moments of Z are 

E[Tr(ZP)] ~ n^-P, Vp ^ 1, 

which is equivalent to say that 

lim K[tr„((nZf )] = 1, Vp ^ 1. 

In all the cases treated above, we leave the proof of the almost sure convergence to the 
reader. □ 

Remark 5.3. Let us notice that for the regimes (/) and (III) studied above, the limit distri- 
bution of the output does not depend on the limit of the input distribution. The result obtained 
in the second regime could have been obtained in a more direct manner, using the powerful 
tools of free probability. For simplicity, let us forget about the normalization of the input 
matrix and notice that the limit distribution of X ®Y is /^(fe), if /U is the limit distribution 
of X and Y is a k x k rank one projector. The partial trace of the randomly rotated input 
matrix is equal to the sum of its k n x n diagonal blocks. Each block is a free compression of 
parameter 1/k (which accounts for a free additive convolution power of k) and the blocks are 
free. Taking the sum of the free blocks explains the other factor k appearing as an exponent 
for the free additive convolution. 



6. Tensor products of quantum channels 

6.1. Motivation and existing results. When studying the question of the additivity of 
minimal output entropies, it is natural to consider products of random quantum channels. 

Before looking in detail at some specific models, let us observe that if one chooses an input 
state which factorizes X12 = Xi ® X2, then 

[$1 (g) ^2]{Xi2) = ^l{Xi) (g) $2(^2), 

and there is no correlation (classical or quantum) between the channels. In order to avoid 
such trivial situations, one has to choose an input state which is entangled. An obvious choice 
(given that ui = n2 = n) is to take X12 = En, the n-dimensional Bell state, and we shall use 
this state in what follows. 

A. Winter and P. Hayden observed that it is relevant in this framework to introduce the 
further symmetry U2 = Ui, as it ensures that at least one eigenvalue is always big. In [9] the 
random matrix inspired by the ideas of Hayden and Winter it was proved that the bounds on 
the eigenvalues could be improved as follows: 

Theorem 6.1. In the regime of k fixed, n 00, the eigenvalues of the matrix Z converge 
almost surely towards: 

• ^ + — p-) with multiplicity one; 

• ~ ■p; with multiplicity k"^ — 1; 

• 0, with multiplicity n'^ — k'^ . 

In the asymptotic regime where n is fixed and k ^ 00, the random matrix Z converges to 
the chaotic state 

P* = —■ 
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If one looks for optimal bounds for the minimum output entropy of $ there is no 
mathematical proof that U2 = Ui is the best choice. Actually, this choice of probability 
measure on U{n) xU(n) does not have full support and one can not rule out that the maximum 
for the minimum output entropy is outside of the support of the probability measure. This 
is what motivates the introduction of the example where Ui and U2 are independent unitary 
matrices. As we will see, this does not yield improvements on the example of Winter with high 
probability. More strikingly, in the regimes that we consider, we will see that the constraint 
U2 = Ui yields no significant improvement to the asymptotic behavior of the von Neumann 
entropies, and this suggests that the simpler random model where Ui and U2 are independent 
could be a candidate for additivity violation with high probability. 

In the forthcoming subsections we analyze both models (independent and conjugate uni- 
taries) in a different asymptotic regime, where both parameters n and k grow to infinity at a 
constant ratio k/n — > c. The model where the quantum channels are independent has received 
less attention from the quantum information community; here, we show that it is intimately 
connected to the (more interesting) case of conjugate channels, by comparing eigenvalue pro- 
files for outputs of channels from the two families. 

6.2. Independent interaction unitaries. Here we consider two independent realizations 
Ui = U and U2 = V of Haar-distributed unitary random matrices on U{nk). For both 
channels the state of the environment is a rank-one projector and we are interested in the 



n 



X random matrix 



where En is the maximal entangled Bell state 



1 " 

-En = - ^ \ei){ej \ (g) \ei){ej 



The diagram associated with the (2,2) tensor Z is drawn in the Figure [71 
We compute the moments E[Tr(ZP)] for all p ^ 1 using the graphical method. We start, as 
depicted in Figure [3 by replacing U* (resp. V*) blocks by U (resp. V) blocks. Notice that 
there are two type of blocks corresponding to the independent random unitary matrices U 
and V (when computing the p-th moment of Z, there are p blocks of each type). This has two 
important consequences: when expanding the diagram in order to compute the expectation 
of the trace, one can only pair U blocks with U blocks and V blocks with V blocks; "cross- 
pairings" between "[/" blocks and 'V blocks are not allowed by the expansion algorithm. In 
addition, one has to index the Weingarten sum by 2 pairs of permutations, one for each type 
of blocks (we shall denote them by au , Pu , cty , (3v G Sp). The 4 permutations are responsible 
for pairing blocks in the following way (1 ^ i ^ p): 

(1) the inputs of the i-th. [/-block are paired with the inputs of the a(/(z)-th U block; 
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(2) the outputs of the i-th f7-block are paired with the outputs of the /3t/(i)-th U block; 

(3) the inputs of the i-th l/-block are paired with the inputs of the av{i)-tli V block; 

(4) the outputs of the i-th F-block are paired with the outputs of the /3y(i)-th V block. 

Since our diagram is made only of unitary matrices (there are no constant non-trivial tensors) , 
the result of the graph expansion is a (sum over a) collection of loops, multiplied by some 
scalar factor. The different contributions of a general quadruple {ajj , [5u , ay , fiv) £ are 
given by (recall that circles correspond to n-dimensional spaces and squares correspond to 
A;-dimensional spaces): 

(1) loops from UU and UU: k*°'u ■ 

(2) loops from oJJ and Uo: n^^'^ 

(3) loops from [/■ and MU : none; 

(4) loops from [/•, •U , V» and n*'-'^u''M- 

(5) loops from DV and VD: A;*"^; 

(6) loops from oV and Vo; n'^^"' "^); 

(7) normalization factors 1/n from the Bell matrices En- n"^; 

(8) Weingarten weights for the [/-matrices: Wg{au(3^^)', 

(9) Weingarten weights for the F-matrices: Wg(ay/3y^). 

Adding all these contributions, we obtain the exact closed-form expression: 

Proposition 6.2. The moments of the random variable Z can be computed as follows: 
(21) E[Tt{ZP)]= A;#"t'+#"^n#(^"'°^)+#(^"'°^)+#(^^''^^)-PWg(ac//3f;i)Wg(av'/3y^)- 

Here we study the asymptotic regime n,k ^ oo, k/n ^ c > 0. Our main theorem is as 
follows: 

Theorem 6.3. Almost surely, the distribution of the output matrix c^n?Z converges towards 
a free Poisson law of parameter . 



Proof. We start by replacing k by cn in equation (j2ip and we obtain 

E[Tr(ZP)] ~ ^ n-'^^c-'^-Moh{au^u^)Mob{avf3v^), 
au,l3u,o!v,l3v&Sp 

where 

Vn = \au\ + \av\ + h~^0!u\ + |7""^av| + \Pu^Pv\ + 2|a[/^^^| + 2|ay/3^^|, 

and 

"Pc = \au\ + + \ctuPu^\ + \oivPv^\- 
Since we are interested in the asymptotic n — > oo (c is a constant), we want to minimize 
Vn- The following inequalities are standard (cf Lemma [23]): 

(22) \au\ + h'^au\;^p-l 

(23) \av\ + h~^av\'^ p-l 

(24) \P^^Pv\,2\auf3u^\,2\av(3v^\ ^ 0, 

and thus P„ ^ 2p — 2 with equality iff au = Pu = OiV = Pv = and a is on a geodesic 
between id and 7. By choosing the obvious rescaling, we get 

1 



lim E 



^ TiUc^n^Z^P^ 



a geodesic a geodesic a(iNC{p) 
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and one recognizes in the last sum the p-th moment of the Free Poisson distribution of param- 
eter c^. This shows that the the matrix Z converges in moments to the hmiting Marchenko- 
Pastur distribution. The argument for the almost sure convergence relies on the Borel-Cantelli 
lemma and can be found in the Appendix. □ 

The von Neumann entropy of the output can be calculated in a fashion simillar to Corollary 
31 



Proposition 6.4. Almost surely, in the limit n 
Z satisfies 

'21ogn- 2^^ + 0(1) 
21og(cn) - ^ + o{l 



oo, the von Neumann entropy of the matrix 



H{Z) 



< c < 1. 



Let us now consider a slightly generalized model of random quantum channels. We intro- 
duce channels $ : — > Aln(C) which have a different tensor product structure at their 
input and output. Here, d is an integer parameter, and we shall always suppose that d \ nk. 
The diagram associated to such a channel is depicted in Figure [H where diamond-shaped la- 
bels correspond to d-dimensional vector spaces and triangle-shaped decorations denote spaces 
of dimension d' = nk/d. The above analysis for a product of independent channels is easily 
adapted to this more general situation: 



E[Tr(ZP)] 



au,l3u,av,Pv'^Sp 



where 



'Pd = IPu^ Pvl and 'Pc = \au\ + \av\ + \auPu^\ + \avPv'^\ 









k. 






■ 


c 


• = 




U 




X 




v* 

• 



Figure 8. A quantum channel with asymmetric input and output tensor structure 



Remark 6.5. If d = d{n) is a function of n such that lim„^oo d{n) = oo, then the considera- 
tions in Theorem \6.3\ carry out to this case and we obtain the exact same limit, a free Poisson 
distribution of parameter (? . The function d = d{n) does not play any role in this situation. 

On the other hand, if the parameter d is constant (inputs of fixed dimension), then the 
limiting behavior changes. Indeed, the minimizing constraint = disappears, and 

the contributing quadruples of permutations become uncoupled: id — > arj = /3(7 — > 7 and 
id — >■ ay = Pv ~^ 7- In conclusion, the asymptotic moments in this case are given by the 
formula, which we summarize in the following proposition 

Proposition 6.6. If d is constant, the limiting distribution of (?r?Z also exists and its limit 
moments are given by: 

-^E[TT((c^n^Z)P)] ~ y c#"^+#°^(i-l"^'"^l. 

id^«j7=/3(7^7 
id— >ay=/3y— >7 
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Question 6.7. We are not able to identify this distribution even though its properties look 
new. We wonder whether this distribution could be related to generalized convolutions of 
Bozejko and coworkers, cf [3] . 

6.3. Conjugate interaction unitaries. To finish, we consider the tensor product of two 
conjugate random quantum channels. As it was emphasized in Section [6. H product channels 
have very interesting eigenvalues statistics and have received a lot of attention in 
the last years because of their usefulness in providing counter examples to different additivity 
conjectures. The purpose of this section is to obtain a description of the behavior of such 
channels in the regime where both n and k grow to infinity at a constant ratio c € (0, oo). 

Hayden and Winter remarked in |19J , that such a conjugate product channel has a very 
important property: the output of the maximally entangled state over the input space has a 
"large" eigenvalue, of size at least l/(cn). The results of [9] show that one expects for this 
model a large eigenvalue Ai = l/(cn) +o(l/n) and (n^ — 1) smaller eigenvalues. The purpose 
of this section is to show that this is indeed the case. Actually, we can prove that the random 
matrix under study has eigenvalues on two scalings: 1/n and In the next theorem, we 

compute the moments of the output matrix Z up to the first order in n. 

Theorem 6.8. Fix some scaling constant c > and consider a sequence of random quantum 
channels ^n,k where n, A; — > oo and k/n — > c. The asymptotic moments of the output matrix 
Z = $ 55 ^{E„) are given by: 

Tr (Z) = 1; 
E TV {{cnZf) = 2 + + 0{n~^); 
E Tr {{cnZ)P) = l + 0{n-^), Vp ^ 3. 

Remark 6.9. Before we prove this result, we would like to draw the attention of the reader 
aware of random matrix theory and matrix integrals, that the symbol 0{n~^) is actually op- 
timal. One can check by inspection that there are terms of order in the expansion of the 
quantities of the theorem. This observation stresses the fact that the matrix model Z does 
not behave like a usual unitarily invariant matrix model, but rather like an orthogonal matrix 
model, even though the underlying invariance group is the unitary group. This technicality 
explains why we can only obtain convergence in probability of the rescaled largest eigenvalue, 
and not the almost sure convergence. 

Proof of Theorem 1 6. 8[ We start from the exact expression at fixed n and k for the moments 
of Z (see [9]): 

(25) E[Tt{ZP)]= A;#"n#("^'')+#(^^)-PWg(a/3-i). 

a,/3eiS2p 

Dropping the number-of-cycles statistics in favor of permutation lengths | • |, replacing 
k ^ cn and using the standard asymptotic expansion for the Weingarten function, we have 

E[Tt{ZP)] ~ c-(l"l+l"^"l)n^'-(l"l+l"^"l+l'^^l+2|"/5-^l)Mob(ar')- 

In order to find the first order asymptotic (in n) of this expression, one has to minimize the 
quantity 

\a\ + |a7-i| + \I35\ +2|a/3-^| 
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over all permutations a, /? € S2p- We start by simplifying this optimization problem over two 
permutations by using the following two inequalities: 

(26) |a| + |a/3-i| ^ 

(27) |a7-i| + |a/3-i| ^ |/37-i|. 

Note that these inequalities can be simultaneously saturated by choosing, for example a = f3. 
So, one is left with the following minimization problem over f3 G S2p'- 

(28) minimize 5i(/?) = + |/?7"^| + 

The main ingredient in tackling this problem is the fact that both permutations 5 and 7 
lie on the geodesic between the identity permutation id and the full-cycle permutation 

7 = (p^...2^1^1^2^...p^). 

This follows from the saturated triangle inequalities \5\ + |(^~^7| = p + p — 1 = 2p — 1 and 
I7l + l7~^7l — 2(1? — 1)-|-1 = 2p — l. If fact, one has 7 = (p'^l^)-^. Under Biane's isomorphism, 
(cf I2.5p . the permutations 6 and 7 correspond to the non-crossing partitions in Figure [9l 



S = 



2T 2B 



1' 1 



B oB 



Figure 9. Non-crossing partitions associated to permutations 5 and 7 

We start by looking at the following simplified minimization problem: 

minimize S2{/3) = \(3\ + \/36-^\ + \/3T'^\. 

Obviously, + |/37~^| ^ 2p — 1, with equality iff /? lies in the geodesic between id and 7. It 
follows from a parity argument that if /? is not an element of the geodesic set id ^ 7, then 
+ 1/37""^ I ^ 2p + 1 and hence, since in this case one has (3^5, S2{(3) ^ 2p + 2. If /? is a 
geodesic element, then S2{(3) ^ 2p — 1, with equality iff /? = 5. 

Since the permutations 7 and 7 are at distance one, the same holds for /?7~^ and /37~^: 

/?7-^=/?7-^-(p^O. 
We have |/37~^| = |/37~^| i 1 and, even more precisely. 




\j3j ""^1 — 1 if p and 1 are in the same block of 
+ l otherwise. 



Note that the condition appearing in the first case can be restated in the following, simpler 
way. It is known that since /3 is a geodesic element, /9~^7 is also on the geodesic, and the non- 
crossing partition associated to /3~^7 is the Kreweras complement of the partition associated 
to /3. The condition that p'^ and belong to the same block of K{P) it is depicted in Figure 
[TOl and it is easily seen to be equivalent to /? ^ 7 (the permutations are compared here via 
their associated partitions). 

It follows that Si{P) = 82(0) a and, in order to conclude, one needs to look at the position 
of the permutation (3 with respect to the geodesic id — > 7. If /? is not an element of id ^ 7, 
then Si{P) ^ 5*2 (/3) — 1^2p-|-2 — l = 2p-|-l which is enough to conclude. We assume from 
now on that /3 is a geodesic element. If |/37~"'^| = |/37~^| + 1, then Si{P) = S2{P) + 1 ^ 2p, 
with equality if and only if /3 = 5 and the conclusion follows. The most difficult case is when 
1/37""^ I = |/37~^| — 1 which is equivalent to the fact that p'^ and 1^ are in different blocks of 
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o* o* o* o* 



Figure 10. Kreweras complement of P 

(3j~'^. We get = 2p — 2 + |/3(5^^|. We claim that for any geodesic (3 such that /3 ^ 7, 

|/?(5^-'^| ^ p. This follows from the fact that the permutation (36~^ = (36 has no fixed points: 
any index is mapped by d to which is then mapped by /3 to some ^ x'^ (the same 
holds for bottom indices). Since it has no fixed points, each cycle of has cardinality 

at least 2 and thus f35~^ has at most p cycles which implies |/3<5~^| = 2p — ^{(36~^) ^ p. 
We conclude that, if a geodesic permutation /3 verifies |/37~^| = |/37~^| — 1, then Si{p) ^ 
2p - 2 + p ^ 2p + 1 for p ^ 3. 

So far, we have shown that the inequality S'i(/3) ^ 2p holds for all /3 and p. Moreover, for 
p ^ 3, we have shown that equality holds if and only if P = 5. For p = 2, using an exhaustive 
search in ^4, one identifies the permutations which saturate the equality: /? G {id, 5, 7}. 

Now that we have completely solved the minimization problem for /3, let us go back to 
equations (j26ti27p and find, for each minimizing P, the values of a which saturate both in- 
equalities. For P = 6, the geodesic id — > 5 has a very simple expression, since (5 is a product 
of transpositions with disjoint support (see the proof of Theorem 6.3 in [9|): 



id^a^(5 44> 3 C A C {1,2, s.t. a = JJ' 



Obviously, one has a6^^ = Ili^Al^"^^^)' ^^'^ thus formula (f27l) reads |a7^"'^| + p — |^| 
Writing explicitly a^~^, one can show that 



#("7" 



1 if yl = 0, 
\A\ otherwise. 



Obviously, A = does not verify the equality, so one is left with 2p — \A\ = \A\ =^ \A\ = p 
and hence a = 6 = p. The other two cases for p = 2 (/? = id and /? = 7) are trivial and yield 
the same result a = p. In conclusion, we obtain for p ^ 3: 

E[Tr(ZP)] = c-Pn-P + o(n-P) 

and for p = 2, 

E[Tr(Z^)] = (1 + 2c-^)n'^ + o{n-^) 
which completes the proof. □ 

At this point, the description of the random matrix Z is not complete: the moment infor- 
mation of the preceding theorem allows us to infer that here are at least some eigenvalues on 
the scale of 1/n and that the rest of the spectrum is distributed on lower scales, such as 1/n^. 
Hayden and Winter's proof of the existence of a large eigenvalue contains, as a byproduct, 
some information on the eigenvector for this particular eigenvalue. Indeed, they use the pro- 
jection on the Bell state to obtain a lower bound for the largest eigenvalue of Z, so one can 
use this projector to obtain more precise information on the eigenvalue distribution of Z. 

In order to obtain information on the rest of the spectrum, we introduce the orthogonal 
projection Q = 1 —E, where E is the maximally entangled state. Using the (rank ii? — 1) 
projector Q, we shall obtain some information on the smallest — 1 eigenvalues of the output 
matrix Z. 
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Theorem 6.10. Almost surely, the matrix (?n?QZQ converges in distribution, to a Marchenko- 
Pastur law of parameter . 

Proof. We compute the moments of the random matrix c^n^QZQ and show that they converge 
to the corresponding moments of the hmit law: 



hm -\ETr{c^n^QZQ)f = [ xPdTT^2{x). 
We start by replacing Q = 1 —E and expanding the product: 

T:t{c^n^QZQ)P = c^^n^P-^E Tr(I -E)Z{1 -E)Z---{1 -E)Z 



= c2Pn2p-2 ^(_i)irHi^)l^-l/-M£^)lETr/(l)Z/(2)Z---/(p)Z, 

where ^ is a set of the 2^ choice functions / : {1, 2, . . . , {I, E}. Notice that in the last 
formula, each Bell projector E is multiplied by a factor — 1/n. 



u 



u 



u 



u 





r ^ 




D 


/ 

/ 


— c 




/ 




D 


/ 

' )C 

V J 


— c 


-5 





u 



3 



Figure 11. Developing Tr(n2QZQ)2 

The moment E Tr f{l)Zf{2)Z ■ ■ ■ f{p)Z is computed with our graphical calculus, and the 
computation is similar to the ones in Theorem I6.3t 

ETr/(l)Z/(2)Z---/(p)Z= ^ A:#"n#("^'"')+#(^^)-P Wg(a/3-i), 

where / G S2p is the permutation associated to the choice function f ^ Tp describing the 
way / connects the different instances of the channel. The exact action of / can be easily 
computed: 



(i-1) 



/ 1(^ + 1)' 



if fii) = I 

if fii) = E, 

if fii + 1) = I 
if /(i + 1) =i^, 



where the arithmetic operations of indices i should be understood modulo p. 

When trying to compute the leading order terms in the expression of ETr(n^QZ(5)P, one 
has to understand the possible cancellations of high powers in n. When writing the exact 
formula for the p-th moment and separating the (a, f3) and the / parts, we get 



-^ETr(cVQZQ)P = c2P ^ 



n 



f<^Fp 
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Note that the sum over f ^ T-p depends only on the permutation a. We show next, that 
for a large class of permutations a (the ones which are responsible for the large eigenvalue 
of size l/(cn) of Theorem 16. Sp . this sum is zero. Let us introduce the set of "vertical line 
permutations" 

V = {cJ G 52p I 3z G {1, . . . ,p} s.t. (j{f) = or a(f) = f} 
= {fj G S2p I (7(5 has at least one fixed point}. 

Fix a permutation a G V and some index i such that a{i'^) = or a(i^) = i'^ . Consider the 
"flip at position i" involution Ti : Tp ^ Tp which maps a choice function / to the function 

I if j = i and /(j) = E, 

Tif -.j ^ (E ifj = i and /(j) = I, 

JU) ifjy^i- 
We shall show that 

(^_i)\r'mn-{\f-HE)\+\af-'\) = (_i)im/)-i{s)i„-{i(T,/)-i(s)i+ia(5^r)"'i)^ 

which will imply that for all a G V, both sums are zero. Since the cardinalities of the sets 
f~^{E) and {Tif)~^{E) differ by exactly one, all one needs to show is that, for all / G J^p, 

\r\E)\ + \ar^\ = mf)-\E)\ + \a{Tjf'\. 
To this end, notice that (the order in which one multiplies the transpositions in not important) 

/= n (a -If J-'') 7 

j ■■ fij)=E 

and hence a{Tif) = af~^ ■ [{i — 1)'^ i^). From this, we find that 

-^77--! I I — 1 if (i — 1)^ and belong to the same orbit of af~^, 

\a[Tij) I = < I I 

I \af 1 + 1 otherwise. 

Let us now suppose that a{i'^) = , the other case a{i^) = being similar. If f{i) = I, 

then f{f) = (i- 1)'^ , af-^{{i - 1)^ = and thus \af-^ \ - \a{T~ff \ = 1. On the other 
hand, f{i) = I {Tif){i) = E and then |/-^(^)| - \{Tify^{E)\ = -1 and one sees that the 
differences compensate. The case f{i) = E is treated in a similar manner. 

We have proved that for all permutations G V, the sum over all choices / G is exactly 
zero; notice that the computations we have done until this point are non- asymptotic, they 
are true at fixed matrix sizes n and k. We now interchange the sums over (a, (3) and /, we 
replace k = cn and we use the first order asymptotic for the Weingarten function: 

irETlic^n^QZQf - V „3p-2-(|/35| + |a|+2|a/3-i|)^2p-(|a| + |a,3-i|) j^^^^. 



^\\f-HE)\„-(\f-HE)\+\af-^\) 

To obtain the dominant power of n, one has to minimize the following quantity over 
(a,/3,/)G(52p\V)x52pX.Fp: 

S{a,l3J) = m + \a\ + 2\al3~^\ + \r\E)\ + \ar\ 
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Since a ^ V, a5 has no fixed point, and hence \a6\ ^ p. Using the facts that |q!/3"^| + \P6\ ^ 
\a6\, \aP~^ \ ^ and \a\ + |a/~^| ^ |/|, we obtain that 

S{a,P,f)^p+\r\E)\ + \fX 

with equality if and only if /3 = a, \a6\ = p and a is on the geodesic between id and /. On 
the other hand, one can easily compute the number of cycles of /: 



#/ 



2 if / = I, 

\f~^{E)\ otherwise. 



Hence, S{a,P,f) ^ 3p — 2, with equality if and only if / = I, /? = a, \a6\ = p and a is a 
permutation on the geodesic id ^ I = 7. Since 7 = 7"^© 7^ is a disjoint union of the two 
cycles 7"^ = {p'^ ■ ■ ■ 2^1^) and 7^ = {1^2^ ■ ■ -p^) = (7"^)~^, the condition that a should be a 
geodesic permutations amounts to a = © , where a^'^ G Sp are geodesic permutations 
with respect to the cycles 7-^'^. One can easily show that i^{{a'^ © a^)d) = #{0^0^), where 
the first permutation is an element of S2p and the second one is an element of Sp. Using this 
equality, the condition |a<5| = p implies a^a^ = idp. When putting all these considerations 
together, one obtains the final formula for the dominant term of the p-th moment of QZQ: 

^ETv{c^n'^QZQ)P ^ c^f-^l^^l Mob(id) = ^ c^*"^. 

id— >a^^7^ id— +cj^— +7^ 

Following the proof of Theorem 16.31 the moments of the Marchenko-Pastur distribution of 
parameter are easily recognized and the converges in moments is settled. The proof of the 
almost sure convergence is more involved and can be found in the Appendix. □ 

From this we deduce the following theorem, which summarizes as the results obtained so 
far in this section 

Theorem 6.11. The eigenvalues Ai ^ • • • ^ A„2 of Z satisfy: 

• In probability, cnXi 1 

• Almost surely, ^■2'_i ^^=2 ^c^n^Xi converges to a free Poisson distribution of parameter 

Proof. Let Ai ^ • • • ^ ^rfi-i be the eigenvalues of QZQ., seen as a matrix in Al„2_x(C). By 
Cauchy's interlacing theorem ([1], Corollary III. 1.5), the eigenvalues of QZQ and those of Z 
are intertwined and satisfy 

Ai ^ Ai ^ A2 ^ • • • ^ A„2_i ^ A„2_i ^ A„2. 



Therefore, the second statement follows immediately from Theorem 16.101 
For the first statement, we have 

1 ^ cnXi ^ c^n^Af ^ c^n^A? + . . . c^n^ A^2 

so the inequality pertains if one takes expectations. In addition, we know from Theorem [67 
that E[c^n^Z3] = 1 + 0(n~^) therefore 

E[cnAi] = 1 + 0{n-^) 

and this proves the first statement. 

□ 

An important result for quantum information theoretic purposes is as follows. 
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Proposition 6.12. Almost surely, in the limit n — > oo, the von Neumann entropy of the 
matrix Z satisfies 



H{Z) 



'21ogn- 2^^ + 0(1) if c^l, 
2 log(cn) - if + o(l) if < c < 1. 

Proof. We use the fact that cnXi ^ 1. Since xlogx ^ — 1 for any x ^ 1, we have 



iXi log(cnAi) ^ (cnXi)^ - 1 ^ ^(cnAj)^ - 1 



en; , „ ^ 

Taking the expectation and using Theorem 16.81 we get 

E[cnAilog(cnAi)] = 0(n"^) 
Similarly, we know by Theorem 16.81 that 

E[Ai] = 0{n-^). 

Putting this together, we obtain 

E[-AilogAi] = o(l). 

We are now left with evaluating 

E[-A2logA2 - . . . - A„2 logA„2]. 

This can be done exactly in the same way as in Corollarv l4.5^ and one then obtains the desired 
formula. □ 

Remark 6.13. Here it is important to remark that the estimate of Propositions [g. 1^ and \6.4\ 
are the same asymptotically. This implies that in this scaling, the choice Ui = U2 is irrelevant 
in the construction of counterexamples to the additivity problem. However it remains to be 
checked whether this scaling indeed yields counterexamples with high probability, and this is 
not clear from our first order asymptotics. 
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Appendix 

In this appendix we present the complete proofs of the almost sure convergence statements 
in Theorems 16.31 and I6.1U[ 

Proof of Theorem 1 6. 31 continued (almost sure convergence). We have already proved the con- 
vergence in moments. To prove the almost sure convergence, it is sufficient to show that for 
each p the series of covariance of the p moments is convergent. A classical application of the 
Borel-Cantelli lemma suffices then to conclude. 

We start with the simplest term, E[tr„2((c^n^Z)P]^. Since we need to compute its first 
two terms in the asymptotic expansion in n, we look at the sub- leading term (n~^) of 
E[tr„2((c2n2Z)P]. Such terms come from permutations for which the exponent Vn has value 
2{p — 1) + 1 = 2p — 1. Analyzing equations (j22ll24p . one sees that the permutations which 
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"almost" saturate the bound are those which verify id ajj = Pu ^ 1^ id 
and |/?[}"^/3\/| = |Q(|}^ay| = 1. In conclusion, we have 



av = f3v ^1 



E[tr„2((c2n2z)Pl 



Taking the square gives 
(29) 



id^a[/=/3y^7 
id— >av'=/3v— >7 



^[tr„2((c2n2z)Pl2 



E 



id^au=Pu=^v=Pv^l 



E 



c 



2#at/ 



id~^au=f3u=civ=l3v^'Y 



E ' 

id^ay=/3y^7 



+ 0{n- 



In order to compute the asymptotic expansion of E (tr„2 ((c^n^Z)^) , one has to consider 

two copies of the diagram corresponding to the p-th power of (?r?Z. The boxes are originally 
connected by the permutation 

72 = (p p - 1 • • • 2 l)(2p 2p - 1 • • • p + 2 p + 1) G 52p. 

After counting the loops, one finds an analogous formula for the mean trace 

,2" 



E 



where 



and 



(tr„,((c2„2z)P) 



E 



n 



4p-4-P„,2„4p-Pc,2 



Mob(a;7/3^^) Mob(Qv/3^^] 



Vn,2 = \av\ + |ay| + |72 + 172 + l^c/Vvl + 2|a;7/3[;^| + 2\ay(iy^\ 



'Pc,2 = |«c/| + \av\ + \oiuPu^\ + loivPv'^l- 
Using the same inequalities and arguments as above, we find that 'Pn.2 ^ 2 [72! 



4p 



with equality iff id ajj = Pu = ctv = /^y — > 72 is a geodesic. Since 72 contains two p-cycles, 
the preceding condition is equivalent to au = /?(/ = ay = Pv = a © a', where a G Sp 
and a' G S{p + l,p + 2, ■ ■ ■ ,2p} ~ Sp are such that id — > a — > 7 and id ^ a' ^ 7 are 
geodesies. Since Mobius functions vanish, this dominating term is equal to the first term 
in the asymptotic expansion (|29|) of E[tr„2((c^n^Z)P]^. The term responsible for the 
contribution comes from permutations au , Pu , ay , Pv G ^2^ such that id au = Pu ^ I2, 



id - 



ay 



Pv 



72 and IPrj'^Pvl 



1. Since, from the geodesic condition. 



a'r, a'lj and ay = a'y a'y, the condition ja^y^'^ayl = 1 is equivalent to either 



a'u 



a'y 



and 



I / // \ — 1 //I 

\{au) ay\ 
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or 

|(a[7)~^ay| = 1 and a'lj = a'y. 
Summing these contributions, one finds the term in from equation (j29p . Hence, the domi- 
nating (n°) and the sub dominating (n~^) terms from E[tr„2((c^n^Z)^]^ and E (tr„2((c^n^Z)P)' 

are equal, which imphes that the general term of the series of covariances has order n^^. The 
series is thus summable and a Borel-Cantelli-type argument finishes the proof of the almost 
sure convergence from Theorem 16.31 □ 

Proof of Theorem \6.1(A, continued (almost sure convergence) . We now prove the almost sure 
convergence statement of Theorem 16.101 We use the same technique as before, showing that 
the covariance series converges. The first step is to analyze the sub- leading terms {n~^) in 
the expression of the p-th moment for one copy of the channel. Recall that the exponent of n 
was given by the expression 

5(a,/3,/) = m + \a\ + 2la/?-'l + \r\E)\ + \ar\ 

Using the triangle inequality |a| + ^ |/|, we split this minimization task into two 
independent problems: 

minimize |/3(^| + 2|a/3~^ |, 

and 

minimize \f~^{E) \ + |/|. 

The 2p—2 minimum in the second problem is reached for / = I; if / is different from I, it follows 
from the above analysis that \f~^(E) \ + \ f\ ^ 2p and thus only / = I contributes to the sub 
leading term. Moreover, a parity argument for the geodesic inequality |a| + ^ |/| 

implies that the permutation a must lie on the geodesic between id and 1 = 7. Let us now 
describe the couples (a,/3) G such that \j35\ + 2|a/?~-^| = p + 1. Since \j35\ + ^ 
\a5\ ^ p, one needs to consider two cases. 

In the first case, we assume that \a5\ = p + \ and a = (3. Since a is a geodesic permutation, 
Q = © and the condition \a5\ = p + 1 is equivalent to = 1. In conclusion, this 

case gives a contribution of 

n ^ 

id— >a"^— ►7"^ 
\aTaB\=l 

In the second case, \a5\ = p and la/?""*^! = 1. This corresponds to = (q^)^^, |a/3^^[ = 1 
and \j35\ = p. Since /3 is at distance one from a, /? = a{i^ f ) for some i, j G {1, . . . ,p} and 
s,t & T, B. Ifs = t, then (3 and thus ^ p which is impossible. We can assume now 
that f3 = a{{^ j^) for some In order to have \P5\ < p, the permutation pS has to have at 
least two fixed points. Using 



and 
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we conclude that, in order to get an contribution, we must have a^(i) = j. Hence, 
for each geodesic permutation a we can find p permutations f3 such that |q;/3~^| = 1 and 
\PS\ = p — 1. We obtain a total contribution of (use Mob(transposition) = —1) 

n ^ 

id— »a'^— ♦7'^ 

Putting the first and the second order contributions together, we obtain the asymptotic 
expansion for the square of the expected normalized trace: 

(30) 



E[tr„2 (c^n^QZQf]^ 



2 

n 



E 

id— * 



E ■ 

id^a^^7^ 



+ 



id^a^— >7''" 



Let us now analyze the second term in the expression of the covariance, E[(tr„2 (c^n'^QZQy)'^] 
The exponent one wants to minimize in this situation is 

S^^\a,(3,f) = + |al +2la/r^| + \f-\E)\ + \af-\ 

where a, /3 are permutations in and the choice function /:{!,... ,p,p+l, . . . , 2p} — > {I, E} 
encodes the way the Z blocks are connected. Note however that in this case, the diagram under 
consideration has at least two connected components, since we are dealing with a product of 
traces. Considerations similar to the ones in the proof of the convergence in moments lead 
to the conclusion that permutations a G V*^^^ do not contribute, so we can restrain our 
minimization problem to the set (^4^ \ V*-^-*) x ^4^ x J^2p- Using the triangular inequality 
|a| + |a/~^| ^ 1/1, we split again our problem in two independent parts: one minimization 
problem for the choice function / and another for the couple (a, P). The minimization problem 
for / is the same as in the single channel case, with the difference that / is defined now on a 
set of cardinality 2p. The quantity |/~^(E')| + |/| is minimized for / = I and the minimum is 
equal to Ap — 4. Notice that in this case, the corresponding permutation I has the following 
cycle structure: 1 = 7 



^T,2 ^B,i ^B,2^ where 
7^'^ = (P^ ip-lf ■■■ 1^) , 



r,2 



r^~ = {{2pf {2p-if 

^B,i = (iB 2^ • . . pB) , 

^B,2 _ ii^ . ^\B 1^ , n\B 



{{p+ir{p+2y 



{2p) 



B\ 



Geodesic permutations id4p share the same cyclic decompositions, and one can easily 

find the dominating term in this case: 

E[{tT„2{c^n^QZQ)Pf] = 



E 



c#""'^+#""'^+o(l), 
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which is the same as the first term in equation ()30p . Let us now move on to the sub- leading 
term in the asymptotic expansion of E[(tr„2(c^n^QZ(5)^)^]. As in the previous case, / = I 
and contributing couples (a, /?) are of two types: permutations such that \a5^'^^ \ = 2p + l and 
a = P or couples such that = 2p — 1 and |a/9~^| = 1. 

The analysis of the first situation is simpler: the cycle structure of the geodesic permutation 
Q implies that |a(5^^^| = 2p+\a^'^a^'^\ + \a^''^a^''^\. Hence, only one of |a^'^a^'^| or |a^'^Q;'^'^| 
is equal to one, the other one being zero. This corresponds to a contribution of (we use the 
symmetry 1 ^ 2 of the problem) 



2 

n 



E 



id^ 
id^ 



^1=1 



The second contribution is calculated in a similar manner to the case of a single trace. 
Permutations /3 at distance one from geodesic a = a^'^ © a^'"^ © a^'^ © a^'"^ such that q^'^ 



{oP"'^) ^ and alpha^'"^ = {oP"''^) ^ are of the form [3 



I. The condition |/3(5(2)| = 2p-l 



implies that we can choose the transposition {i^ j^) and that [a^'^ © a'^'^](i) = j. This last 
condition implies that i and j have to be in the same half of the set {1, . . . + 1, . . . , 2p} 
and, using again the symmetry between the first and the second trace, we can write the final 
contribution: 



n 



id- 



E 



id- 



Summing the leading {n^) and the sub-leading (n ^) contributions and comparing to equation 
(130]) . we find that 

,-2\ 



E[(tr„2 (c^n^gZQf)^] -E[tr„2 (c^n^QZQf]^ = 0{n-'^), 
and the convergence of the covariance series follows, ending the proof. 



□ 
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